Entanglement of Fock-space expansion and covariance in light-front Hamiltonian 
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We investigate in a model with scalar "nucleons" and mesons the contributions of higher Fock states 
that are neglected in the ladder approximation of the Lippmann-Schwinger equation. This leads to 
a breaking of covariance, both in light-front and in instant-form Hamiltonian dynamics. The lowest 
Fock sector neglected has two mesons in the intermediate state and corresponds to the stretched 
box. First we show in a simplified example that the contributions of higher Fock states are much 
smaller on the light-front than in instant-form dynamics. Then we show for a scattering amplitude 
above threshold that the stretched boxes are small, however, necessary to retain covariance. For 
an off energy-shell amplitude covariance is not necessarily maintained and this is confirmed by our 
calculations. Again, the stretched boxes are found to be small. 



I. INTRODUCTION 

In his famous article of 1949, Dirac described a num- 
ber of ways how to set up a framework of Hamiltonian 
dynamics. Two of these are most important. In instant- 
form (IF) Hamiltonian dynamics one specifies the initial 
conditions on the equal-time plane t = 0. Of the ten 
Poincare generators six are kinematic, i.e., do not con- 
tain interactions, and are therefore conserved quantities, 
and four are dynamical: the three boost operators and 
the energy. In light-front (LF) Hamiltonian dynamics 
one uses the coordinates 



A' 



(A^,^^), A^ 



A°±A' 



(1) 



and quantizes on the plane ~ 0. Then only three 
operators are dynamical. Two of them involve rotations 
around the x or y-axis. Therefore LF time-ordered am- 
plitudes are not invariant under such rotations. 

The question of rotations in LF dynamics (LFD) was 
discussed before with the aim of constructing the angular 
momentum operators, see e.g. Fuda |^], and the review 
by Carbonell et al. While these authors emphasize 
the algebraic properties of the generators of the Poincare 
group, we stress in the present paper the connection be- 
tween expansions in Fock-space and covariance. It has 
been remarked before by Brodsky et al. that the higher 
components in Fock space contribute to the difference 
between the Bethe-Salpeter equation and the evolution 
equation in LFD. These authors do not give numerical 



*VUTH 98-20, submitted to Phys.Rev.C 



^nico@nat.vu.nl 

*blgbkkr@nat. vu.nl 

^ karmanov@sci.lebedev.ru 



estimates of the corrections. The latter has been done 
by Mangin-Brinet and Carbonell |^ , and by Frederico et 
al. who studied the same model and found the ef- 
fect of higher Fock states on the binding energy to be 
small. In a calculation of positronium, Trittmann and 
Pauli used an effective theory, where the effects of all 
Fock states are included in the interaction. They found 
rotational symmetry to be restored in the solution. 

In this paper we consider first standard LF quanti- 
zation and discuss the problem of noncovariance, which 
includes violation of rotational invariance, in the frame- 
work of LF time-ordered perturbation theory. We give 
numerical results for the simplified case of two heavy 
scalars exchanging light scalar particles. This choice is 
motivated by the popular meson-exchange models in nu- 
clear physics. We do not include the internal spin degrees 
of freedom, as this is a complication that may obscure the 
main point of our investigation: the connection between 
the breaking of covariance and a truncation of the expan- 
sion in Fock space. In two interesting papers, Fuda 
reported on detailed calculations of realistic one-meson 
exchange models in both LF and IF dynamics. There the 
emphasis is on the comparison between the two, when in 
both cases the ladder approximation is made. It is the 
purpose of the present paper to show to what extent the 
ladder approximation violates covariance. 



A. Suppression of higher Fock states 

A reason why LFD is preferred by many is that higher 
Fock states are said to be more strongly suppressed in this 
form of dynamics. As a disadvantage the lack of man- 
ifest rotational invariance, and therefore covariance, is 
mentioned. We call a symmetry manifest when it is con- 
nected to a kincmatical operator. Then all time-ordered 
diagrams exhibit this symmetry. Equal-time ordered di- 
agrams lack boost invariance whereas on the LF the lon- 
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gitudinal boost P+ is a kinematical operator. Therefore, 
if one complains about lack of manifest covariance, one 
should include not only rotational invariance but also 
other nonmanifest symmetries. One reason why people 
have rather stressed rotational invariance comes easily to 
mind: in many cases it is easy to convince oneself by 
inspection whether a matrix element is rotational invari- 
ant, viz when the amplitude can be expressed in terms 
of scalar products of three-vectors. On the other hand, 
it is not more difficult to test numerically for invariance 
under boost transformations than for invariance under 
rotations. Indeed, the method used in the present paper 
can easily be extended to check for boost invariance. 

A way to test for covariance is to compare the LF 
time-ordered diagrams to the covariant amplitude, since 
we know that the latter is invariant under any of the 
Poincare symmetry operations. For on energy-shell am- 
plitudes (S'-matrix elements) there is an exact equality, 
as was proved by Ligterink and Bakker and which 
is confirmed in our results. Off' energy-shell there is a 
deviation, which, however, is found in the present paper 
to be surprisingly small. 

So why are we using these LF time-ordered diagram 
in the first place, when there is an equivalent covariant 
method available? We do, because we want to determine 
the properties of the bound state using the Hamiltonian 
form of dynamics. In this method, covariance can never 
be fully maintained. However, one may try to apply it 
in such a way that breaking of covariance is minimal. 
In many applications in nuclear physics a one-meson ex- 
change approximation is made for the interaction and 
the scattering amplitude is computed by formally iterat- 
ing this interaction, leading to the Lippmann-Schwinger 
equation in the ladder approximation. In this approxi- 
mation one retains two- and three-particle intermediate 
state and neglects Fock states containing four or more 
particles. These Fock sectors are needed to make the 
sum of LF time-ordered diagrams equal to the covariant 
amplitude, exhibiting the symmetries under all Poincare 
transformations. If these contributions are large, one can 
expect a significant breaking of covariance, since the LF 
time-ordered diagrams are only invariant under applica- 
tion of the kinematical symmetries. 

For this reason we concentrate in this article on the 
determination of the contributions of these higher Fock 
states. Our main concern shall be the box diagram. Then 
we label the correction as 71^+ . We shall calculate 71^+ 
explicitly for the box diagram with scalar particles of 
different masses. The box diagram can be associated 
with the two-meson exchange between two nucleons. If 
spin would be included several well-known complications 
would arise, the most important one being the occur- 
rence of instantaneous propagators |jl^,|ll|. We do not 
want these complications to interfere with the main point 
of our investigation: the connection between Fock-space 
truncations and lack of covariance. Therefore spin is 
omitted. We have not included crossed box diagrams, 
because they are not relevant for a discussion on covari- 



ance, since both the crossed and noncrossed box diagram 
are covariant by themselves. 



B. Setup 



First we explain the Lippmann-Schwinger formalism 



III 



and the special role of the box diagram. In Sec 
we describe how to calculate both the covariant and the 
LF time-ordered amplitudes. After this, we are ready 
for our numerical experiments. In Sec. |y| the masses of 
the external particles are chosen in such a way that on- 
shell singularities of the intermediate states are avoided, 
and therefore it is easy to compare IF and LF Hamilto- 
nian dynamics. In that section it is shown that 7?.4+ is 
much smaller in LFD than in IF dynamics (IFD), con- 
firming the claim that in LFD higher Fock states are more 
strongly suppressed. Moreover, it tells us that covariance 
is more vulnerable in IFD than on the light-front! 

After this exercise, we concentrate on the LF, and in 
Sec. V] we calculate the LF time-ordered diagrams for 



the more interesting case in which we have particles of 
fixed masses m (called nucleons) and fj, (called mesons). 
As the process we are concerned with, scattering, is above 
threshold, we have to deal with on-shell singularities. We 
show that the breaking of covariance is again small. 
Although in Sec. VII, where we discuss off-shell ampli- 



tudes below threshold, no on-shell singularities are en- 
countered, matters become more complicated because 
the notion of the cm. frame becomes ambiguous, since 
the total momentum is dynamical and found to be un- 
equal to the combined momentum of the two particles, 

+ q^. However, we are still able to relate the breaking 
of covariance and Fock-space truncation. 

The lack covariance of the LF time-ordered amplitudes 
means that the amplitude depends not only on the scalar 
products of the external momenta, but on the angles be- 
tween the quantization axis and the external momenta 
as well. Consequently, the amplitudes must have singu- 
larities as a function of these angles in addition to the 
familiar singularities as functions of the invariants. The 
positions of these singularities are found analytically in 
Sec. VIH, in the framework of explicitly covariant LFD. 



This gives a qualitative understanding of the numerical 
results in Sees. ^ and VI. In Sec. IX explicitly cova riant 
LFD is applied to the off energy-shell results of Sec. VH . 



II. THE LIPPMANN-SCHWINGER FORMALISM 

The Hamiltonian method aims at the determination 
of stationary states, i.e. eigenstates of the Hamiltonian. 
Here we take the Yukawa-type model with scalar coupling 



5$2 



(2) 



Two types of particles are considered: "nucleons" (N, <&) 
with mass m and "mesons" (m, (f>) with mass fj.. The 
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Hamiltonian H = P~ consists of a part Hq which de- 
scribes free particles and a part V which describes the 
interaction: 



H = Hn 



V. 



(3) 



We shah denote the second term on the r.h.s. as the po- 
tential. The problem of constructing the Hamiltonian 
from the underlying Lagrangian has been recently re- 
viewed by Brodsky et al. [1^] . Here we study two-nucleon 
states only. Moreover, we neglect self-energy diagrams. 

We consider i/o, the kinematic part of the Hamilto- 
nian in the two-nuclcon (2N) sector. In the instant-form, 
quantization is carried out on planes of constant time 
{equal-time planes). Then we find for two particles of 
mass TO and momenta p and q respectively 



(4) 



which leads to both negative and positive energy solu- 
tions. It is well-known M] that in this form the overall 
momenta and the relative momenta are difficult to sepa- 
rate. 

In LF quantization the square root, and therefore the 
negative energy solutions are absent. The interaction- 
free part of the two-body Hamiltonian is 



Hq - — 



-|- TO 



2p^ 



2q^ 



(5) 



which demonstrates that positive energies occur for pos- 
itive plus-momenta. Moreover, one can easily separate 
the motion of a many-particle system as a whole from 
the internal motion of its constituents in the LF case 

We shall focus on light-front quantization of our model 
in which the interaction of the nucleons is due to meson 
exchange. We write the potential in the form 



V ^Vi + V2 + V3 + 



(6) 



where the subscript denotes the number of mesons simul- 
taneously exchanged. The potentials only contain irre- 
ducible diagrams to prevent double counting. Vi contains 
one-meson exchanges only: 



V, = 



(7) 



The irreducible diagrams contributing to Vi are depicted 
in Eq. (|^). In these diagrams time goes from left to right. 
The nucleons are denoted by solid lines, and the mesons 
by dashed lines. Irreducible diagrams contributing to 
V2 are those diagrams of order that cannot be sepa- 
rated into two pieces by cutting two nucleon lines or two 
nucleon lines and one meson line only. In terms of Fock- 
space sectors this means that Vi contains two-nucleon 



and one-meson intermediate states, and V2 contains only 
two-nucleon two-meson intermediate states. 

The potential Vi is a covariant object in the case the 
external lines are on shell. The meaning of the equal 
sign in Eq. (j^) is that the full covariant amplitude can 
be written as a sum of two LF time-ordered diagrams. 
Whereas the Feynman diagram contains the propagator 
l/((q' — q)"^ — the LF time-ordered diagrams con- 
tain the energy denominator 1/{P~ — Hq), P~ being the 
parametric energy. Hq is the sum of the kinetic energies 
of the particles in the intermediate state: 



Ho 



E 



TO, 



2fc+ 



(8) 



The two diagrams contain ^^-functions of the plus- 
component of the momentum of the exchanged meson: 
one has the factor 9{p~^ — q'^), the other 0{q~^ —p'^)- 

In a Feynman diagram the external lines are on mass 
shell and the initial and final states have the same energy, 
which coincides with the parametric energy. Then the 
minus component of the total four-momentum of a two- 
particle state satisfies the relation 

P =p +q =p +q = \ (9) 



2p^ 



2q^ 



As the minus component of the total momentum is the 
only dynamical momentum operator, the other three 
components are conserved in any LF time-ordered dia- 
gram. For instance, P"*" = p~^ + q'^ = p'~^ + q'^ . This 
conservation law is very important in LF quantization. It 
leads immediately to the spectrum condition: in any in- 
termediate state all massive particles have plus-momenta 
greater than zero and the sum of the plus-components of 
the momenta of the particles in that state is equal to the 
total plus-momentum. 

The expansion in Fock space does not coincide with an 
expansion in powers of the coupling constant. This can 
easily be seen when one considers an approach closely 
resembling the Lippmann-Schwinger method. The eigen- 
states IV' > of the Hamiltonian 



H\iIj >= p-\tP >, 



(10) 



are also solutions of the Lippmann-Schwinger equation 
1 



P--Hq 

The formal solution of this equation is 



1=0 



P- -Hn 



V] \6> 



(11) 



(12) 



An equation similar to Eq. (|ll| ) exists for the scattering 
amplitude: 



V 



V 



(13) 
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If one substitutes Vi for V in these equations one obtains 
the ladder approximation. This approximation does not 
generate all diagrams, so one needs to add corrections. 
At order this correction is V2 : 



(14) 



If one takes into account all the contributions to V from 
Eq. (pi) then the full scattering amplitude is 



boost invariance is not manifest. Therefore, we refer to 
breaking of covariance, which is a general feature of any 
form of Hamiltonian dynamics, if one truncates the the 
Fock-space expansion. 

We would like to estimate the contribution of the miss- 
ing diagrams, irrespective of the strength of the coupling. 
It is not possible to do this in a completely general way, 
so we perform our numerical calculations for the box dia- 
gram only. We assume that our results will be indicative 
for the higher orders too. 

We define the fraction 



(15) 



In the ladder approximation one only takes Vi into ac- 
count. Effectively, one then describes the full interaction 
between two nucleons by 



(16) 



In this approximation intermediate states containing 
more than three particles do not occur. This implies 
that time-ordered box diagrams with four particles in 
the intermediate state are neglected, as we can see if we 
compare the expansions in Eqs. ( p^ and (p^). As the 
individual diagrams contributing to V2 are not covariant, 
the sum of box diagrams produced by the ladder approx- 
imation is not covariant. 

Using equal-time quantization, twenty out of the 
twenty-four possible time-orderings have intermediate 
states with more than three particles. On the LP, the 
spectrum condition kills many of the time-ordered dia- 
grams. There are six nonvanishing diagrams, of which 
four only contain two and three-particle intermediate 
states. One concludes that the one-meson exchange ker- 
nel neglects the majority of the contributing time-ordered 
box diagrams in equal-time quantization, whereas on the 
LP most of the nonvanishing diagrams are taken into 
account. This does not mean necessarily that in IP dy- 
namics the ladder approximation misses most of the am- 
plitude, since the missing diagrams have smaller magni- 
tudes. The contribution of the missing diagrams needs 
to be investigated in order to see how much the higher 
Pock sectors are suppressed. 

There is one thing which seems to complicate mat- 
ters on the LP. The individual LP time-ordered diagrams 
are not rotational invariant. When a number of them is 
missing the full amplitude will also lack rotational invari- 
ance, as is mentioned often in the literature. This fea- 
ture does not occur on the equal-time plane, since here 
rotational invariance is a manifest symmetry. However, 
in other types of Hamiltonian dynamics other symme- 
tries are nonmanifest. In IP Hamiltonian dynamics, e.g.. 



(17) 



The subscript 4 indicates that this variable includes all 
diagrams having at least 4 particles in some intermedi- 
ate state. Por TZ^^ one would have to add the diagrams 
containg five and six-particle intermediate states in the 
numerator, as these give nonvanishing contributions in 
the instant-form. The diagram in the denominator is the 
covariant diagram. 

We shall show that the correction V2 is indeed much 
less important numerically in LFD than in IPD. So the 
2N2m-state is in LPD much less important than in IPD. 
We conjecture that this property of LPD that the Pock- 
state expansion converges much more rapidly than in IPD 
persists in higher orders in the coupling constant. 



III. THE BOX DIAGRAM 

In the previous section we saw that the lowest level 
at which breaking of covariance is to be expected is the 
two-meson exchange diagram, also referred to as the box 
diagram. Implicitly, in Sec. [| we took the particles to be 
scalars, as we did not include instantaneous terms which 
are related to spin- 1/2 particles. Although a bound state 
of scalar particles is not found in nature, we do not in- 
clude spin because we want to avoid in this investigation 
the complications due to instantaneous terms. 

We look at the process of two nucleons with momenta p 
and q respectively, coming in and exchanging two meson 
of mass /i. The outgoing nucleons have momenta p' and 
q'. The kinematics is given in Fig. g. 



FIG. 1. Kinematics for the box diagram. The arrows de- 
note the momentum flow. 

The internal momenta are 
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fci = k, 

fc2 = fc - p', 

h ^ k - p - q, 
ki — k ^ p. 



(18) 
(19) 
(20) 
(21) 



We have to keep in mind that these relations only hold 
for those components of the momenta that are conserved. 



A. The covariant box diagram 



The covariant box diagram is given by 

-i d^'fc 



(22) 



Mi: 



where the subscript "Min" denotes that the integration 
has a Minkowskian measure. The imaginary parts ie of 
the masses have been dropped. 



B. Equivalence 

If the external states are on energy-shell, that is, 

P- =p- +q- =p'- +q'-, (23) 

then the time-ordered diagrams are the same as those 
derived by integrating the covariant diagram over light- 
front energy k~ . In that case we have 



(24) 



The example of the box diagram with scalar particles of 
equal mass has been worked out before by Ligterink and 
Bakker M. 



C. The LF time-ordered diagrams 

It is well-known [|ll],0 how to construct the LF time- 
ordered diagrams. They are expressed in terms of inte- 
grals over energy denominators and phase-space factors. 
In the case of the box diagram we need the ingredients 
given below. The phase space factor is 



(25) 



Without loss of generality we can take p+ > The 
internal particles are on mass-shell, however, the inter- 
mediate states are off energy-shell. There are a number 
of possibilities for the intermediate states. We label the 
corresponding kinetic energies according to which of the 



internal particles, labeled by fci . . .k^ in Fig. |l|, are in 
this state. 



Hii = q- 



k^ 

2fc+ 



Ik A 



H13 — 



ki + k^ + 



2kt 



H12 = q'- 



kf + 



2fc+ 



Hzi =p - 



H2A = q' +p 



2fc^ 



2/c^ 



^23 = P 



fej + Kg 



2k^ 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



A minus sign occurs if the particle goes in the direction 
opposite to the direction defined in Fig. |l[ All particles 
are on mass-shell, including the external ones: 



q = 



P 



q-^ + 
p + m 



q = 



q^ + m 



P 



2q'+ 

/ I 2 9 

p^ + 



(32) 



2p+ ' ' 2p'+ 

We can now construct the LF time-ordered diagrams. 
The diagrams (33) and ( ]35| ) will be later referred to as 
trapezium diagrams, (|^ as the diamond and (^) as the 
stretched box. 

^ ' /o {P-- Hu) {P-- H^z) (P-- i?i2) 

/ = / d^fc 



-27r dfc+ 



d^fc^ 



d^fc^ 



p,+ 4> {P-- Hu) {P-- His) {P-- H23) 

-27r dfc+ 

+ (P-- Hs4) {P-- H13) {P-~ H23) 

-2n dk+ 

'p'+ {P-- Hu) (P-- H2i) {P-- H23) " 
= 0. (37) 



(33) 
, (34) 
, (35) 
, (36) 



The factor 2t: matches the conventional factor in 
Eq. (|2^). The last two diagrams are zero because we 
have taken p^ > p'^ and therefore these diagrams have 
an empty fc+-rangc. If wc take p"*" < p'"*", which case will 
also occur in upcoming sections, the diagrams (^) have 
nonvanishing contributions. 



IV. A NUMERICAL EXPERIMENT 

We look at the scattering of two particles over an angle 
of tt/2. In Fig. |^ the process is viewed in two different 
ways. 
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z 





(a) 

FIG. 2. (a) Two particles come in along the a;-axis. They 
scatter into the y — z plane over an angle of 7r/2. The az- 
imuthal angle is given by a. (b) Another viewpoint. The 
outgoing particles go out along the y-axis. One can vary 
the orientation of the light-front vector uj with respect to the 



Fig. ||a pictures the situation where the scattering 
plane is rotated around the x-axis. The viewpoint in 
Fig. concentrates on the influence of the orientation 
of the quantization plane, and is connected to ex plicitly 
covariant LFD, as will be discussed in Sec. VIII. Both 



as will be discussed in Sec. VIII . 
viewpoints should render identical results, since all angles 
between the five relevant directions (the quantization axis 
and the four external particles) are the same. We choose 
for the momenta 









0, 


0), 


(38) 






~v\ 


0, 


0), 


(39) 






0, 


— w^. 


-«^), 


(40) 


e = 




0, 






(41) 



indicating that we have chosen the fixed quantization 
plane — (Fig. ^). The incoming and outgoing 
particles have the same absolute values of the momenta 
in the cm. system. Therefore 



= („-)2 = („!^)2 + („.)2 

The Mandelstam variables are 



s 

t 

u 



ip- 



„'^2 



- -m\ 



(42) 



(43) 
(44) 
(45) 



We are now ready to perform the numerical expe riments 
for three cases, which are described in Sec. |V- VII. In the 
experiments two parameters are focused on. We shall 
vary the azimuthal angle a in the y-z-plane, 



a = arctan — , 
and the incoming c.m.s. -momentum 



(46) 



V. LIGHT-FRONT DYNAMICS VERSUS 
INSTANT-FORM DYNAMICS 

One of the claims of LFD is that higher Fock states are 
more strongly suppressed than in IFD. We can investigate 
this claim for the box diagram easily in the following case. 

We take the external states on energy-shell (p3|), such 
that the equality ( p^ holds. At the same time we avoid 
on-shell singularities for the intermediate states by giving 
the external particles a slightly smaller mass m', 



2 2 

p < m , 



(48) 



such that we can still relate the amplitude to a 5-matrix 
element. 

The process we look at is described in the previous sec- 
tion and has two scalars of mass m' coming in along the 
X-axis, interacting, and scattered over a scattering angle 
of 7r/2. Stretched boxes give maximal contributions (see 
next section) if the quantization axis is in the scattering 
plane, which is the case if the azimuthal angle a = tt/2. 
TZ^+ is easily found by calculating the stretched box. 



TZ^^ +TZl^, however, has twenty nonzero 



contributions. As an example, we show the six contribu- 
tions to TZ^ in Fig. ||. 



FIG. 3. Time-ordered diagrams that contribute to TZs. 
The diagrams in the first column have five particles in the first 
intermediate state. The diagrams in the second column have 
five particles in the last intermediate state, and the diagrams 
on the right have five-particle intermediate states for both the 
first and the third intermediate state. 

This helps us to understand why TZ^^ = 0. All con- 
tributing diagrams contain vacuum creation or annihila- 
tion vertices, which are forbidden by the spectrum con- 
dition. 

There are twelve diagrams contributing to TZq, and all 
contain vacuum creation or annihilation vertices. There- 
fore TZq^ vanishes. 

We calculated TZ^^ by subtracting the four diagrams 
only containing three particle intermediate states from 
the full sum. This sum can be obtained by doing the 
covariant calculation, or by adding all LF time-ordered 
boxes. Our results are given in Fig. H. We also calculated 
TZ^^ 



v = v''. (47) 

In our experiments we will omit the units for the masses, 
which are MeV/c^. 
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Sec. VII] we identify the singularity /sing, 
the nonstretched boxes as 





FIG 



100 200 300 V 400 

4. LF time-ordered boxes for a scattering angle of 

7r/2 as a function of the incoming momentum v. We also give 

the ratios of boxes with at least four particles (7^4+ and 7^4+ ) 

or five particles (7?.r+ , TZ^f — 0) in one of the intermediate 

states. 



We conclude that on the LF contributions of higher 
Fock states are significantly smaller than in IFD. In the 
limit w — > the ratio TZ^^ goes to zero, because the phase 
space becomes empty. However, in IFD there is a finite 
contribution of TZ^^ = 4.5% in this limit. Even if one 
includes five-particle intermediate states the LF is the 
winner by far. 

Note that m', given by Eq. (48), varies as a function 
of p^, and therefore also as a function of w, but is inde- 
pendent of a. The deviation of m' from m is small; less 
than 2.3% for v < 200, and less than 9% for v < 400. 
As the deviation of the mass m' from m is only small, 
we are convinced that these results are indicative for cal- 
culations above threshold. However, we do not want to 
do these calculations, because then one needs to subtract 
the on-shell singularities of the equal-time ordered boxes. 



VI. NUMERICAL RESULTS ABOVE 
THRESHOLD 

As in Sec. we look at the scattering of two particles 
over an angle of 7r/2. We focus on LFD, and therefore we 
simply write TZi — Tl^+ ■ We do no try to avoid on-shell 
singularities by taking different masses for the internal 
and external nucleons. Two nucleons of mass m — 940 
scatter via the exchange of scalar mesons of mass fi = 
140. Again, the vertex is a scalar and no spin is included. 



A. Evaluation method 

Contrary to the case considered in Sec. now it 
is not straightforward to evaluate the contributions of 
the LF time-ordered boxes, because the nonstretched 
boxes c ontain on-shell singularities, thoroughly analyzed 
in Sec. VIII . Here we briefly sketch how we deal nu- 
merically with the singularities. Using the analysis of 



(Tk I 



d3fc(/ - /si„g) + / d^fc L 



and rewrite 



(49) 



The integrand /sing has a simple algebraic form, such 
that the integration in one dimension over the singular- 
ity can be done analytically, and the remaining integral is 
regular. This integral is then done numerically by math- 
EMATICA. The integral over (/ — /sing) was implemented 
in FORTRAN. These two numb ers are ad ded to give the 
results presented is Sees. VI B and VI C^ . 



B. Results as a function of a 

We shall now vary the direction of v' , given by the az- 
imuthal angle a, however, not its length. Therefore the 
Mandelstam variables are independent of a, and we ex- 
pect the full amplitude to be invariant. We tested this 
numerically for a number of values of v. In the region 
< a < TT we used the formulas ( p3|) until (Bq) . In the 
region tt < a < 27r the diagrams (p4|) and ( pq ) vanish. 
However, then there are contributions from the diagrams 
in (|37|). The results are shown in Fig. ||. The results are 
normalized to the value of the covariant amplitude. The 
contributions from the different diagrams vary strongly 
with the angle a. Since the imaginary parts are always 
positive, they are necessarily in the range [0, 1] when di- 
vided by the imaginary part of the covariant amplitude. 
The real parts can behave much more excentric, espe- 
cially for higher values of the incoming c.m.s. -momentum 
V. An analysis of the a-dependence is given in Sec. VIII. 



Clearly the LF time-ordered diagrams add up to the co- 
variant amplitude, so we see that in all cases we obtain 
covariant (in particular rotationally invariant) results for 
both the real and the imaginary part. 



C. Numerical results as function of v 

We look at scattering in the a;-2;-plane (a — 7r/2), be- 
cause in that case the contributions from the stretched 
boxes are maximized. The results are shown in Fig. ^. 

We see that the stretched box, the diagram with two 
simultaneously exchanged mesons, is relatively small at 
low energies, but becomes rather important at higher en- 
ergies. We depict the ratio of the stretched box to both 
the real part and to the total amplitude. Since the real 
part has a zero near v = 280, the first ratio becomes 
infinite at that value of the incoming momentum. 
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(a) Real part for v = 40: TZ^ = 0.07%. (b) Imaginary part for v = 40. 




7r/2 TT 3n/2 a 27r 

FIG. 5. Amplitudes above threshold from a = to a = 2tt. TZa is the maximal fraction of the stretched box to the absolute 
value of the sum. 
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FIG. 6. Real (a) and imaginary part (c) of the LF time-ordered boxes above threshold for a = n/2 as a function of the 
momentum of the incoming particles v. The inset (b) shows the ratio of the stretched box to the real part of the amplitude 
(7^4^) and to the absolute value (TZi). 



VII. NUMERICAL RESULTS OFF 
ENERGY-SHELL 

In the previous section, we tested covariance of the LF 
formalism for amplitudes with on energy-shell external 
particles, by using the c.m. frame, where P-^ = p'-^ + 
q'^ = and — p'^ -\- q'^ — 0. However, on the LF the 
operator is dynamical, and the last equality does not 
hold anymore off energy-shell, as one can easily verify in 
the following way. Consider the case of a bound state 
with mass Ai < 2m. The bound state is off energy-shell, 
and therefore we have 

p-<p'-+q'-. (50) 

The plus and transverse momenta are kinematic, so 

P+=p'+ + q'+, (51) 



P^ 



p'^ + q'^. 



(52) 



Adding Eqs. (|50|) and (|5l]) gives 
P' < p'' + q" 



(53) 

q"' > 0. 

Therefore the two outgoing particles cannot have exactly 
opposite momenta as in Eqs. (HO) and (|4l|). In terms 



If = 0, then Eq. (53) imphes that p' 



of the explicitly covariant LED, introduced in Sec. VIII 



this reflects the fact that the off energy-shell relation be- 
tween p' + q' and P contains extra four-momentum like 
in Eq. (p9|). What was the reason that we chose opposite 
momenta in the previous sections in the first place? Our 
reason was that we wanted to have a manifest symmetry 
of the amplitude, because it is obvious that the Man- 
delstam variables s, t and u remain the same under the 
rotations we investigated. We investigate how the off- 
shell amplitude breaks covariance by using momenta for 
the outgoing particles, such that the Mandelstam vari- 
ables remain constant, at the same time satisfying the 
conditions (|ol)-(|3l). 

As in the previous sections, we shall fix the direction 
of the incoming particles, as in Eqs. ( ^8| ) and (^), and 
vary the direction of the outgoing particles. Therefore, 
by construction, s is invariant. In order to guarantee the 
invariance of the other Mandelstam variables t and u we 
must take p-q' and p-p' constant, even off energy-shell. 
These inner products are: 

(54) 



p-q' =p^q' 



p -q 



p^ 

nl 



2(1- Xp,)P^ 



+ il-Xp,)P^ 



p-p =pp 



p'+p 



p-p 



P^ 



,P^ 



p -q 



(55) 



p-p 
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We have introduced the fractions 



(a) V = 40: 7^4 = 0.11% and TZoS = 0.08%. 



(56) 



Since the perpendicular momenta are conserved we have 
in the cm. system p'^ — — g'^, so the inner products of 
the perpendicular momenta can be written as 



P^-p!^= \p^\\p'^ 
pW^ 



cos 9, 
b^lb'^lcos^, 



(57) 
(58) 



where 9 is the scattering angle. We can now solve 
Eqs. (|§]5^) for Xp> , \p'^\ and 9. There are many 
curves satisfying these conditions. For uniqueness we de- 
mand that the curve goes through the point in which 
Xp, = Xp = 1/2, = \p^\ and 9 = n/2. We find that 
the curve is then parametrized by 



9 = Tr/2, 
7— r = 4(p +m ). 



(59) 
(60) 



From Xpi and p''^ we determine p'"*" and p'~ , and finally 
also p'^. Because the particles come in along the x-axis, 
the above relations define an ellipse in the y-z-plane. In 
the case of IFD these ellipses would reduce to circles with 
center of the origin and radius v. 




100 

1 




\ -100 

'\ 

\\ -100 


100 i 


(b) 





Possible momenta p' and 



tides (fat line) in the scattering plane (horizontal: p'^ 



of the outgoing par- 
and 

q''" , vertical: p'^ and q'^) for two cases: (a) w = 40 and (b) 
V — 200. The momentum p' + q' is indicated by the dot. As a 
reference we have drawn the locus for on energy-shell external 
particles (thin line): a circle centered at the origin. 

In Fig. 1^ we have indicated the y and z-components of 
the momenta of the external particles for the two cases 
we investigate. We see that Eqs. (^) and (|5^) hold. The 
off energy-shell momenta form an ellipse. The deviation 
from a circle with radius v is hardly visible. 

In Fig. [s] we show the contributions of the different 
boxes and their sum as we vary the angle a. The calcu- 
lations are the same as in the previous section, using the 
formulas (|3^)-(|3^), except that the momenta of the out- 
going particles have changed such that ( pof ) and (^^ are 
satisfied. As there does not exist a covariant amplitude in 
the off energy-shell case, we normalized the curves shown 
by dividing the amplitudes by their sum at a = 7r/4. 



0.8 



0.6 




3n/4 a ^ 2n 



(b) V = 200: 7^4 = 2.04% and 7^off = 5.64%. 
sum 




7r/4 7r/2 3n/4 2k 

FIG. 8. The LF time-ordered boxes as a function of a. 

In Fig. ^ we see the consequences of the off energy-shell 
initial and final states. Condition ( |2^ ) is violated, there- 
fore Eq. (|^ does not hold and a breaking of covariance 
can be expected. We see that the contributions of the 
higher Fock states TZi are smaller than the effect of the 
off-shellness TZqS ■ This is confirmed in Fig. ||, in which a 
is fixed and the incoming c.m.s-momentum v is varied. 




200 







100 200 300 V 400 

FIG. 9. The LF time-ordered boxes as a function of v for 
a = 7r/2. The inset shows the maximum contributions of the 
stretched box, and the maximal breaking of covariance. 
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From Fig. || we infer that the full amplitude is maximal 
at a = and a = it. The minimum is reached at a = 
7r/2 and a = 37r/2. Therefore the maximal breaking of 
covariance of the amplitude can be calculated by taking 
the difference of the total sum at the angles a = and 
a — 7r/2. We see that at typical values for incoming 
momentum [v < 50) TZoS is small, even smaller than 
TZ^. However, at higher momenta it dominates over the 
stretched box. In this region we see that the stretched 
box contributions remain small. 

A detailed explanation of the results for off energy-shell 



amplitudes is given in Sec. IX 



VIII. 



ANALYSIS OF THE ON ENERGY-SHELL 
RESULTS 



The angle dependence of the LF time-ordered ampli- 
tudes found numerically, can be understood analytically. 
The variation of the LF amplitudes with the angle a 
means that they have singularities in this variable, either 
at finite values of a or at infinity. They should disappear 
when they are summed to give the covariant amplitude. 
These singularities can be most conveniently analyzed in 
the explicitly covariant version of LFD (see for a review 
1^). In this version the orientation of the light-front is 
given by the invariant equation lu-x — 0. The amplitudes 
are calculated by the rules of the graph technique ex- 
plained in Ref. S. After a transformation of variables, 
these amplitudes coincide with those given by ordinary 
LFD. However, they are parametrized in a different way. 
The dependence of the amplitudes on the angle a means, 
in the covariant version, that they depend on the four- 
vector CO determining the orientation of the LF plane: 
M = M{p, q,p', q' , uj) Hence, besides the usual Man- 
delstam variables s ( p3| ) and t (Q) the amplitude AI de- 
pends on the scalar products of lu with the four-momenta. 
Since to determines the direction only (the theory is in- 
variant relative to the substitution lu — > auj), an ampli- 
tude should depend on the ratios of the scalar products 
of the four-momenta with ut. Hence [n5 : 



M — M{s, t, Xp, Xp'), 



where 



u>-p 



u-{p + q) ' 



iCp' 



UJ-p 



uj-{p + q) 



(61) 



(62) 



The formulas (|6^) coincide with the definitions ( pq ) if 
we use the z-axis as the quantization axis. The lu- 
dependence is reduced to two scalar variables Xp and 
Xpi , since the direction of lu is determined by two an- 
gles. Hence, this amplitude should have singularities in 
the variables Xp and Xp'. Their positions will be found 
below. The amplitude corresponding to the sum of all 
time-ordered diagrams should not depend on Xp and Xpi . 

Let us find the physical domain of the variables Xp 
and Xp' , corresponding to all possible directions of lu for 



fixed s,t. In the cm. system the variables Eq. 
represented as: 



uj-p, 



uj-p 



3) are 



(63) 



where, e.g., ili-p is the scalar product of the unit vectors 
LU — Lu/\uj\ and p = p/\p\ in three-dimensional Euclid- 
ian space, and v = s/A — rrfi is the momentum of the 
particle in the cm. system. The Eqs. ( ^3| ) determine 
an ellipse in the Xp-a;p'-plane. Its boundary is obtained 
when LU is in the scattering plane, that is n-p = cos 7 or 
fi-p = cos(7 — 6), where 7 is the angle between p and 
LU in coplanar kinematics and 6 is the scattering angle in 
the c.m.s. The case when Q is out of the scattering plane 
corresponds to the interior of the ellipse. For a scatter- 
ing angle = 7r/2 the ellipse turns into a circle, shown in 
Fig. 




-"-max -"-p' 



FIG. 10. Physical region in the a;p-Xp; -plane. 

For the kinematics shown in Fig. ^ and Eqs. (|3|)-(|lj), 
when w _L p, it follows from Eq. (^) that the value 
'p — ^, whereas for a given a we obtain: 



i.e 



Xp is fixed: x 



2vo 



sma. 



(64) 



with vq — \J m? -\- v^, varies along a straight line when d) 
is rotated in the y-z-plane. The bounds of the physical 
region of Xpi are: 



V 

2^' 



When < a < tt/2, Xp' moves from 1/2 to x^nin- When 
TT /2 < a < n, Xp' moves in the opposite direction in the 
same interval. This explains why all the curves in Figs. ^ 
and ^ in the interval < a < tt are symmetric relative 
to a = tt/2. 

When n < a < 3tt/2, Xp' moves from 1/2 to Xmax and, 
finally, when 37r/2 < a < 2tt, it goes back in the same 
interval. As in the previous paragraph, this explains why 
all the curves in Figs. || and ^ in the interval n < a <2tt 
are symmetric relative to a — Stt/2. When a = 'k/2 
and 37r/2, the values of Xp' are on the boundary of the 
physical region. 
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Note also that the amphtudes for the trapezium 
(dashed and sohd curves in Figs. ^ and |^) are evidently 
obtained by the replacement p q, p' q' , which, 
according to the definition in Eq. (|62|), corresponds to 
Xp' — > 1 ~ Xpi . This is the same as the replacement 
a ^ 2tt — a in Eq. (|64|). Therefore the curves for the 
other trapezium, when a goes from 27r to 0, repeat the 
curves for the trapezium, when a increases from to 2tt. 
The same is true for the other diagrams (diamonds and 
stretched boxes). 



A. Trapezium 

The method of finding the singularities of the LF dia- 
grams was developed in ||l5|| . Here we restrict ourselves 
to the example of the diagram (^3|). Its counterpart in 
the explicitly covariant LFD is shown in Fig. 



COX2 



COX:! 



COX, 



p k p 

FIG. 11. The trapezium in explicitly covariant LFD. 

The dotted lines in this figure are associated with fic- 
titous particles (spurions), with four- momenta propor- 
tional to Lu. The four- momenta of the particles (the spuri- 
ons not included) are not conserved in the vertices. Con- 
servation is restored by taking into account the spurion 
four-momentum. In the ordinary LF approach this cor- 
responds to nonconservation of the minus-components of 
the particles. The spurions make up for the difference. 

According to the rules of the graph techniques jsj , one 
should associate with a particle line with four-momentum 
p and mass m the factor 9{uj-p)5{p'^ — w?) and asso- 
ciate with the spurion line with four-momentum wr^ the 
factor 1/(t,; — ie). Then one integrates, with measure 
A'^ki/ {2tiY , over all the four-momenta ki not restricted 
by the conservation laws in the vertices, and over all dr^. 
The expression for the amplitude of Fig. ^ is: 



[uj-k)6{k^ - m^) 

x0{ujip' - k)) d{(p' -k + ujTif - n^) 
x9{lo-{p + q-k)) 6{{p + q - k + ujt^Y - m?) 
x9{lu-{p - k)) 5{{p ~k + uoT-if - 
dri dT2 dT3 d*/c 



(ti - ie) (t2 - ie) (T3 - ie) (27r)3 



(66) 



Like in Eq. (p2|), we omit the coupling constant. Perform- 
ing the integrals over and dfco in Eq. ( |66| ) by means of 
the (5-functions, we get: 



e{ujip'-k)) 9{uj-{p + q - k)) 



0iuHp^k)l 



d^fc 



2£fe(2^)3- 



(67) 



By transformations of variables (see for the details ap- 
pendix B of [|j) the expression Eq. (|67| ) can be trans- 
formed such that it exactly coincides with Eq. (^3|). 

For Feynman amplitudes the method to find their sin- 
gularities was developed by Landau in ||l6|| . A method 
very similar to that one can be applied to time-ordered 
amplitudes. If we would omit for a moment the 9- 
functions in Eq. ( p7| ) and would not take into account 
that k^ = m^, for finding the singularities we should 
construct the function (pi — ai(/i^ — {p' — fc)^) +a2{m'^ — 
{p + q — fc)^) -I- 03(777,^ — {p — fc)^) formed from the de- 
nominator of (|67|). The singularities of the trapezium 
are found by putting to zero the derivatives of (pi with 
respect to ai_3 and with respect to k. However, the 
trapezium may have singularities corresponding to a co- 
incidence of the singularities of its integrand with the 
boundary of the integration domain caused by the pres- 
ence of the ^-functions. So, we must find a conditional 
extremum. The restrictions can be taken into account us- 
ing Lagrange multipliers [p5| . Hence, we should consider 
the function: 

if = ai{i? - {p' - kf)+a2{n? - {p + q~kf) 
+ a:i{^? -{p- kf) + ai{e - m") + 7iw(fc - p'), (68) 

where 04 and 71 are the Lagrange multipliers. One 
should also consider the functions obtained from Lp at 
Ci\ = 0, then at ai — 0, then a-i = 0, then at ai = a2 = 
0, etc. One should not consider the function obtained 
from (|6^) by = 0, since the integral (|67| ) contains the 
three-dimensional integration volume d^^ and the con- 



dition k^ 



cannot be deleted. Therefore there is 



no need to introduce the term 72^-^, since the fc^ = m? 
condition prevents this term from being zero and, hence, 
does not impose any restrictions. The case 71 = repro- 
duces the singularities of the Feynman graph. Therefore 
below we shall consider the case 71 7^ resulting in the 
singularities in the variables Xp,Xpi. We suppose that 
uj-p > Lo-p' . This corresponds to the condition p+ > p'"*" 
of Sec. IIL In the kinematics shown in Fig. |2|this means 
that Xp' < 1/2 and < a < tt. In this case, the second 
and the third 0-functions in Eq. ( |67| ) do not give any re- 
strictions and can be omitted. Therefore we omit also 
the term j3U!-{p — k) + ^^uj-{p + q — k). 

The derivatives of ip with respect to fc, the a's and 71 
give: 



dip/dk ^ ai2{p' - fc) + a22{p + q - k) 
+ a^2(j) -k)+ ai2k -f 71a; = 0, 



with 



{p'~kf^^?, (p + q-kf^ 



{p~kf 



j-k — LU-p' . 



(69) 



(70) 
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We multiply Eq. (|69|) in turn by (p' — fc), {p + q — k), et cetera, and get: 



(69) xp' ~k 
(6£) y. p + q- 
(6£) X p — fc : 
(6£) X fc : 
(6£) X u : 



ai/i^ + a2'irn? + a^^^ + a^{s — 2m^) + 7i(l — Xp/) — 0, 
ai{2y? — t) + a2fJ,^ + a^^^"^ — a^^i'^ + 71 (xp — Xp') = 0, 
—ai^^ + 0:2(5 — 2m?) — a^^J? + 2Q!4m^ + 71 Xp' 
02(1 — Xp') + a3(xp — Xp') + a4Xp' ~ 0. 



0, 



(71) 



These equations have a nontrivial solution iff: 



2//2 




(2^2 _ 









2to2 


(s- 2m2) 


(l-Xp') 






2/^2 




{Xp Xp 


(s - 2m2) 




2rr? 


Xp' 





(1 - Xp') 




Xp' 






(72) 



Eq. d?^ ) is quadratic in Xp' . Its solution is simple but 
lengthy. We show it for the particular case of the kine- 
matics of Fig. H supposing that the particles in the cm. 
system have momenta v. In this case s and t are given 
by Eqs. (ph and (E|). The solution of Eq. (H) 



is: 



V 



1 

x", — - ± 

^ 2 2vq + 8/^2^2 + 4^4 



(73) 



The position of singularity in the variable Xpi is denoted 
by Xp'. According to |16| the behavior in the vicinity of 



should be either logarithmic, proportional to |xp 
Xp, or proportional to |xp' — x^, |^ log(xp' — x^,), where 
/3 is a noninteger number. 

When 11 V we find from Eq. (|7^): 



1 V , u . 
2±2^(^ + £^)- 



(74) 



Comparing with Eq. (|65|), we see that at small fi or at 
large u the singularities come closer to the physical region 
of Xp' . We will see below that this will be a property of 
all the singularities depending on fi and v. This explains 
the numerical results, showing that with an increase of 
V the graphs of the amplitudes versus a become more 
sharply peaked. 

Now consider the case when one of the a's is zero. Let 
ai = 0. Then Eq. (|6^) is reduced to: 



(fi = a2{nn? — {p + q ~ k)"^ 
+ a^{k^ 



asi^i^ ~ {p ~ kf) 
n^)+jiuj-{k-p'). 



(75) 



Similarly to the previous case we get an equation, which 
can be obtained from Eq. (^2|) by deleting the first row 
and column: 



2m2 
(1 



2to2) 

Xp' ) 



A^^ 
2fi^ 

(Xp - 



s — 2m 
2 



2m2 



(1- 

{Xp - 
Xp' 





= 0. 



(76) 



Under the given kincmatical conditions its solution with 
respect to Xp' is 



1 _l_ /X ^ 4^2 



(77) 



In the limit fj, —^ or v 00 these singularities are again 
approaching the physical region. 

Let a2 = 0- The equation is obtained from Eq. ( |7^ ) 
by deleting the second row and column: 



2^2 {2n^-t) -^2 

(2^2 _ 2/^2 _^2 (-^^ _ 



-A* 



2to2 Xb' 



('^p "^p' ) '^p' 

Its solution reads: 







= 0. 



(78) 



X 



Q Xp^(4m2^ ~ fi^ — fit± 2\/t^tm^ + — 4to2/x2) 



In the limit 



4m2^2 - (i - ^2)2 

00 it is simplified: 
,2 



X, 



4z;2 2^2 ' 



(79) 



(80) 



Let as = 0. The equation is obtained from Eq. ( [T^ ) 
by deleting the third row and column: 



2^2 ^2 _^2 

^2 2m2 (s-2to2) (1-Xp') 

—fi"^ (s — 2rn2) 2777,2 

(1 — Xp') Xp' 



= 0. 



The determinant in Eq. IKW) can be evaluated: 



4sXp' — 4sxp' 



4m2 - ^72 = 0. 



The solutions of Eq. (82) are: 

1 . ^i;2 + /Lt2/4 



± 



2vo 



(81) 



(82) 



(83) 



For /i ^ they also approach the boundary of the phys- 
ical region of Xp' . 
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Now consider the cases when a few coefficients are zero. 
Let ai = as = 0. The equation can be obtained from 
Eq. ^72) by deleting the first and third rows and columns: 



2m^ 



(s — 2m ) 2m^ 

(1 -Xp/) 



This equation reduces to: 



(s-2m2) (l-Xp/) 







= 0. 



Ep/ — Xp'S + TT? = 0. 



Its sohitions are: 

Q _ 1 V 
' " 2 + 2^ 



1 V 
2^ 2v^ 



(84) 



(85) 



(86) 



Since we consider the interval < a < tt corresponding 
to Xmin < Xpi < 1/2, the singularity at x^, = Xmax is 
beyond the physical region, whereas the singularity at 
Xp, = Xniin is just On the boundary of the physical region. 
The amplitude in this point gets an imaginary part: 



Im 



Im 



^ at Xpf ^ Xi^ 



at Xp' — x^i 



(87) 



Eq. ( pq ) corresponds to a = 7r/2. This explains, why 
all the dashed curves of the imaginary parts in Fig. || go 
through zero at the point a = 7r/2. 

Now put ui — a2 = 0. The corresponding equation is 
obtained from Eq. ( [Z^ ) by deleting the first and second 
rows and columns: 



2/x^ — (xp 

—^^ 2iv? Xpi 

{xp Xp' ) Xp' 



Xpi) 



Eq. ( p9| ) reads: 

2 2 /'r, 2 2\ I 2 2 

Xp,m — Xp'Xp(2m — ^ ) + XpUi 
Its solution is: 



(89) 



(90) 



(91) 



These two singularities are fixed points in the complex 
plane. At Xp = 1/2 and /i <C m they are approaching the 
point Xpi — 1/2 in the physical region, i.e., a = and 
a = TT. 

The case a2 = as = leads to the equation obtained 
from Eq. ( fz^ ) by deleting the second and third rows and 
columns: 



2/^2 

— /i^ 2m^ Xp' 
Xp> 



0. 



(92) 



It gives Xp, = 0. This is a fixed singularity in the non- 
physical region. 



Above we have considered the region cu-p > wp' . In 
the region uj-p < u-p' the integration domain is restricted 
by the step function 9{iLj-{p — k)) instead of 9{lu-{p' — k)) 
in Eq. (p?)). The integrals defining these amplitudes de- 
fine different analytic functions depending on the region 
considered. In the point Xpi — 1/2, i.e. at a = and 
a — IT, the values of the functions coincide, but their 
analytic behavior is different. 

This can indeed be seen in Fig. |[ The slopes at a = 
and a = TT are different. 



B. Diamond 

The diamond corresponding to Eq. (Q) is shown in 
Fig. 111. 



P k P 

FIG. 12. The diamond in explicitly covariant LFD. 

The analytical expression is: 



= / 9{uj-k)6{k^ - m2) 



X6'(w(fc - p'))S{{k -p' + LJTl - UJT2f - M^) 

x9{uj-{p + q~ k))S{{p + q - k + ^72^ - -rr?) 
x9{ujip - k))6{{p -k + cuT^f - fi^) 
dTidT2dT3 d'*fc 
(ri - ie){T2 - ie){T3 - ie) {2tt)^ 



(93) 



Performing the integrations in Eq. ( |93| ) over r^, we get: 



9iLoik-p')) 



^2 _ (fc _ p,)2 + _igl^_Pl_[^2 _ (p + 5 _ k)^] 



9{uj\p + q^k)) 9{LU-{p-k)) 
•m? — {p + q — kY — (p — A;)2 



9{uj-k) 



2£fc(27r)3 



(94) 



We still suppose that to-p > uj-p' . However, now, in 
contrast to the trapezium, uj-p' < uj-k < uj-p, and both 
restrictions have to be taken into account. 

In order to find the singularities, one should consider 
the extremum of the function: 
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[ uj-{p + q-k) 

+a2 {ni^ - (p + q- kf} +03 {^J■^ - {p - kf} + ai{k^ - w?) 

+-/iujik -p') +-/2Loik-p). (95) 

At Lo-p' = uj-p, i.e. at a = and a — n the integration domain vanishes and the diamond becomes zero, as shown in 
Fig. 1^. It remains zero in the interval tt < a < 27r. 

C. Stretched box 



The stretched box, corresponding to Eq. (36) is shown in Fig. |l 



- COT2 



P "m^ P 

FIG. 13. The stretched box in explicitly covariant LFD. 



The analytical expression is: 



[LO-k)S{k^ - m'^) 

xe{uj-{k - p'))5{{k -p' + 0JT2 - wrs))^ - m") 
x9{uj-{p + q~ k))S{{p + q - k + UTi + tJTa — 1^X2)^ — m^) 
xe{uj-{p - k))5{{p -k + uT^f - 
dTidT2dr3 d'^fc 



(ti - ie)(T2 - ie)(r3 - it) (27r)3 
Performing the integrations in Eq. (Bq) over r^, we get: 



e{Loip + q-k)) 



{p + q-kf + ""''J^llp,^^ if^' -{k-p'r] 



,^-{k- p'Y + ^^j^^W ~{P- kf] " " '''^^'"^ 

In order to find the singularities, one must consider the extremum of the function: 

Lp = a-Y\m ~{p + q~k) H - [k - p ) ) 

uj-{k-p) 

2 ^i\2 , ^-(k-p') 2 2^ 



(96) 



+a2 - (fc - p'Y + - b - kY) 

+a■i{^Y - (p - kf} + ai{k^ - m^) + -injj-{k - p') + 72t^-(fc - p). (98) 

Calculating the derivative of Eqs. ([zsl), (|9^ ) and (|8|), for example, with respect to ai-4, at 71 = 72 = 0, one finds 
identical equations determining the singularities. These do not depend on Xp,Xpi and coincide with the ones of the 
Feynman graph. Similarly one can see that any singularity depending on Xp and Xpi cannot appear in a separate 
diagram only. It appears at least in two amplitudes. These singularities cancel each other in the sum of the amplitudes. 
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IX. ANALYSIS OF THE OFF ENERGY-SHELL 
RESULTS 

The off energy-shell amplitude is shown graphically in 
Fig. 



COT 



cox' 



(a) (b) (c) 

FIG. 14. Off energy-shell amplitudes in explicitly covariant 
LFD: (a) The trapezium, (b) The diamond, (c) The stretched 
box. The external momenta are the same for all diagrams. 

It contains incoming and outgoing spurion lines with 
the momenta lot and wr', respectively. The conservation 
law has the form: 



P- 



- LOT 



P 



P. 



(99) 



From Eq. (|99| ) one can infer that if P = 0, then p' +q' ^ 0, 
as it was indicated for the z-components in Sec. VI] . To 



parametrize the off energy-shell amplitude, we introduce 
different initial and final Mandelstam variables s: 



l\2 



and the total mass squared: 



M" = {p + q — ojt) = [p + q — ujt') 



l\2 



(100) 



(101) 



So, in general, the off energy-shell amplitude is 
parametrized as: 



M = M{s,s' ,M'^,t,Xp,Xp>). 



(102) 



The on e nergy-shell amplitude Eq. ( |61| ) is obtained from 
Eq. (102) by the substitution s — s' — M'^. One can 
also consider the half off energy-shell amplitude with one 
incoming or outgoing spurion line. It is obtained from 
Eq. (102) by the substitutions s = ^ s' or s' = 

M'^ ^ s. 

In the case of the trapezium, Fig. the external spu- 
rion lines enter and exit from the diagram at the same 
points as the momenta p and p' . So, they can be incor- 
porated by the replacement: 



p p ~ LOT, p 



P 



(103) 



This corresponds to new masses of the initial and final 
particles for the bottom line: 



mf — {p — lot)'^ = — Xp{s — A^^), 
i} = {p' -ujT'f =n? -Xp\s' - JV^). 



(104) 



With t hese ne w masses, one can repeat the calculations 
of Sec. VIII A and find the singularities of the off energy- 
shell amplitude for the trapezium. The masses of the 
intermediate particles are not changed. 



For other diagrams, both for the diamond and the 
stretched box, in contrast to the trapezium, the spurion 
line enters in the point where the momentum g' go out 
from the graph. This means that the calculation has to 
be done with the following external mass of this particle: 



(105) 



whereas the mass of the particle with momentum p' is m. 

Like in the case when all masses are equal, the sum 
of all time-ordered graphs with masses different from 
the internal ones, but the same in all the time-ordered 
graphs, would not depend on lo. However, now we take 
the sum of the graphs with different external masses in 
different particular graphs. This sum cannot be obtained 
by the time-ordering of a given Feynman graph. In this 
case the w-dependence is not eliminated in the sum of 
all the graphs, and the exact off energy-shell amplitude 
in a given order still depends on lo. An example of this 
dependence is shown in Fig. |[ 

The off-shell amplitude is not a directly observable 
quantity. It may enter as part of a bigger diagram. 
Therefore, the off shell amplitude may depend on lo. This 
w-dependence is not forbidden by covariance and, hence, 
does not violate it. On energy-shell this dependence dis- 
appears. 



X. CONCLUSIONS 

If sufficient caution is exercised, invariance of 5'-matrix 
elements can be maintained in Hamiltonian formulations 
of field theory. A necessary condition to be fulfilled is 
that all Fock sectors included in the Feynman diagrams 
that contribute to a perturbative approximation of the S- 
matrix are retained. For the specific case of scalar field 
theory at fourth order in the coupling constant we have 
determined the magnitude of the breaking of covariance if 
only the diagrams generated by the ladder approximation 
to Hamiltonian dynamics are included. The remaining 
terms, the stretched boxes, were found to contribute a 
small fraction, less than 2% for small to intermediate 
c.m.s. momenta, of the total amplitude. This fraction is, 
however, increasing with energy. 

It was found, in a calculation closely approximating the 
first one, that the violation of invariance is much larger 
in instant-form dynamics than in light-front dynamics, 
confirming quantitatively what has been claimed in the 
literature. 

In both cases we determined quantitatively the depen- 
dence of the six LF time-ordered diagrams on the orien- 
tation of the light-front. We verified that, although the 
individual diagrams depend strongly on the orientation, 
their sum does not, as it should not. This dependence of 
individual diagrams may be interpreted as a breaking of 
rotational invariance. 
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Having established numerically that invariance of the 
S'-matrix elements is obtained only if all Fock sectors 
relevant to a certain order in perturbation theory are 
included, we extended our investigation to amplitudes 
that are off-energy-shell. Such amplitudes are not S'- 
matrix elements, calculated between asymptotic states, 
from — oo to -|-c» in time. They are elements of an S'- 
matrix calculated for finite light-front time, i.e., defined 
on a light-front in the interaction region, not moved to 
±00 Therefore they depend on the orientation of this 
light- front. They either occur as parts of larger diagrams 
that are invariant, or in the calculation of LF wave func- 
tions. Not being invariant, the sum of the six LF time- 
ordered diagrams corresponding to the box is expected to 
depend on the orientation of the light-front. We found 
the variation of the sum of these six diagrams to grow 
more strongly with increasing relative momentum than 
the fraction carried by the stretched boxes. 

All these results point to the conclusion that for low 
and intermediate momenta, e.g. those relevant for the 
bulk of the deuteron wave function, the higher Fock com- 
ponents are very small and are expected to play a minor 
role in LF dynamics. We conjecture that this conclusion 
remains essentially valid for higher orders in perturbation 
theory. 

Two remarks are in order here. First, if bound or scat- 
tering states at high values of the relative momentum 
are to be calculated, the higher Fock states will become 
much more important. Secondly, in the present work we 
neglected spin. It remains to be seen to what extent the 
special effects of spin, notably instantaneous propagators, 
will influence our conclusions. 

A final point concerns the dependence of the individual 
diagrams on the orientation of the light- front. By an 
analysis very close to the Landau method for Feynman 
diagrams, we were able to explain all the peculiarities of 
the angular dependence in terms of the occurrence and 
position of singularities of the time-ordered diagrams as 
a function of the angles and their locations. In particular 
the symmetries of the angular dependence and the cusps 
showing up at specific orientations could be explained 
fully. 
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